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Abstract 

We calculate the branching ratio of S ^ J/^' vr" with a mixed formalism that combines the 
QCD-improved factorization and the perturbative QCD approaches. The result is consistent with 
experimental data. The quite small penguin contribution \n B ^ J /'^ tt^ decay can be calculated 
with this method. We suggest two methods to extract the weak phase (3. One is through the 
dependence of the mixing induced CP asymmetry Sjj^i^o on the weak phase/3 , the other is from 
the relation of the total asymmetry Acp with the weak phase /?. Our result shows that the 
deviation ASj^^,^o of the mixing induced CP asymmetry from Sin{—2j3) is of 0(10"^) and has 
much less uncertainty. The above 0(10"^) deviation can provide a good reference for identifying 
new physics. 
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B physics is entering the era of precision measurement, It is not far from reveahng new 
physics beyond the Standard Model(SM). Many authors have studied the topics and suggest 
some windows for looking for new physics(NP)[l|-[9|. Because falvour- changing neutral 
current (FCNC) processes only occur at the loop-level in the SM , so they are particularly 
sensitive to NP interactions. It was pointed out that 5° — mixing and decays are good 
places for new physics to enter through the exchange of new particles in the box diagrams, 



or through new contributions at the tree level 



m- 



12| , so 5° — system has been studied 



in many papers for probing new physics 13|, ll7|. B J/^tt*^ decay is a good mode for 



looking for new physics and extracting the weak phase (3 . The direct CP asymmetry Cj/^t^o 
and the deviation AS'j/^^o = Sj/^^^o — sin(— 2/3) of the mixing -induced CP asymmetry 
from sin(— 2/5) in this decay arise from quite small penguin contribution in the SM, so these 
quantities are sensitive to new physics effect. Comparing the prediction of CP asymmetry 
in the SM with the experimental data, one can find new physics signal. Thus it is essential 
to calculate the ASj/^t^o and Cj/^j^o in B ^ J /'^ti^ in the SM accurately. 

The deviation ASj/^^^o = Sj/^^o — sin(— 2/3) or direct CP asymmetry Cj/^t,o in i? — >• 
J/\E' 71^ decay have been studied in Ref. ISj by fitting to the current experimental data, the 
result is Cj/^t^o = 0.09 ± 0.19 which has very large uncertainty. In that case we can not say 
anything about new physics effects. 

In order to reveal new physics effects, we need both better theoretical prediction and 
experimental measurement with less uncertainties. That is the aim of our present paper. 

In what follows, we first evaluate the penguin pollution effect by_a method which have 
been used to explain many B decays into charmonia successfully 19|, l20|] . We find the penguin 
pollution in the B J rc^ decay is quite small, the deviation ASj/^t,o = S'j/^^o— sin(— 2/5) 
in B ^ J/^ 7r° decay is O{10~^), which means that the measured deviation ASj/^t^o at 1% 
will indicate the presence of new physics. 



The latest experimental data of ASj/^t^o is Sj/^t^o = —0.4 ±0.4 2l|, which has large error, 
so we are expecting to have more precise measurement in the near future. 
The decay rate of of i? — ^ J/\I' ir^ can be written as 

^ GUl-rl + U-rlMr (1) 



327rmB '2 

with r2 = mj/^/niB, n = m^/niB. 

The amplitude A consists of factorizable part and nonfactorizable part. It can be written 
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FIG. 1: Nonfactorizable contribution to the BS' — > J jil) vr*^ decay 



as 



A = A 



NF 



A 



VERT 



A 



HS 



(2) 



where Anf denote the factorizable contribution in Naive Factorization Assumption(NF), 
Avert is the vertex corrections from Fig. [Il(a)-(d) , Ahs is the spectator correction from 
Fig.[Il(e)-(f). 

The factorizable part Anf in Eq- ([2]) for i? ^ J/'^tt^ decay can not be calculated reliably 
in the pQCD approach, because its characteristic scale is around 1 GeV. We parameterize 
the sum of the factorizable part Anf and the vertex corrections Avert as, 



Anf + Avert = aeffmlfj/^Ff^''{m'^j/^){l - rj) 



(3) 



where fj/^ is decay constant of J/ip meson. 

For the B — > it transition form factors, we employ the models derived from the light-cone 
sum rules 



22( 1 . which have been parameterized as 



+ 



r2 



1 — q'^/rrii 1 — q'^/m 



(4) 



fit 



with ri = 0.744, r2 = —0.486, mi = 5.32Gev, Tn^fn = 40.73G'ef for S — vr transition. 



The factorization and vertex correction from Fig. [Il(a)-(d) can be calculated in the 
QCDF 23|. Summing up the factorizable part and vertex correction , we can get the Wilson 
coefficient Cg//, 



^^ + ^^iv: + 4^iv:^^ 



.18 + 121n^ + /, 



-v: 



C3 + Tf + 



Cr a. C 



Ca -18 + 12 In 



nib 



rrib 



fi 



+0^ + ^ + -r^Ce 6 - 121n— - /, + C7 + r^f + a 



CiO 



Nr in N, 



Nr 



Nr 



with the function, 



fi 



2y^r 



3{1-2X2) , „ 2^ , 2ri(l-X2) 

In X2 — 3m + 3 ln(l — H ^ 

1 — X2 1 — r2X2 



(5) 



(6) 



The spectator corrections Ahs from Fig. [TJ(e)-(f), can be calculated reliably in the pQCD 



as in Ref . 
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2Q|, 



(7) 



where the amplitudes A^^"'/'^'^'* and J^^^^'^''^^ result from the {y — A){y — A) and {y — 
A)(y ^ A) operators in the effective Hamiltonian, respectively. Their factorization formulas 
are given by the pQCD approach. In the calculation of A^i'^/'^^^ and m!^^^'^'^^ ^ because J/V' 
is heavy, we reserve the power terms of r2 up to C(r|), the power terms of ra up to . 

16TTmlCF\f2Nr [dx] bidbi<^B{xi , h) 

x{ [(1 - 2rl + rl){l - X2)$.(X3)^^(X2) + ^(r^ - 4)^^{x3)<il\x2) 
-r,(l - rl)x3<^P^{x3)'^ l{x2) + {^rUl - X2) + (1 - r^Xs) $^(x3)*^(x2) 

X Ei^4 (4^^ ) /li^^ (Xi , X2 , X3 , 61 ) 

- {x2 - X2r\ + X3 - 2r2X3 + rlx2,)x2,)^T,{x-i)'^^ {X2) 



+Tl{2r^^l{x^) - -(1 - r2)<l>^(x3))^*(x2) 

-r,(l - r2)x3$^(x3)^'L(x2) - (2r2x2 + (1 - rl)x^ <l>t(x3)^^^(x2) 



xEi,4(tJ^)/lJ^(Xi,X2,X3,6i) , 

lQTcmlCF^2Nc [dx] bidbi^B{xi,bi] 



(8) 



X { [(1 - X2 + rjx2 + X3 - 2r^x3 + r^Xg - r^)$^(x3)^'^(x2) + 
r2(2r^$t(a:3) - ^(1 - rl)<^^{x,))^\x2) 

-r^(l - r2)x3$^(x3)^^(x2) - (2r2(l - x^) + (1 - r2)x3) $t(a;3)^^(x2) 
X EG{t^j'^)h^j\xi, X2, xs,bi) 

;i _ + r2^)x2$.(x3)*^(x2) + ^(r^ - r2^)r2<l'^(x3)**(x2) 



-r^(l - r2)x3$^(x3)^^(x2) + (2r2x2 + (1 - rl)x3) <l>t(x3)*^(x2) 
xE6(tf (^i,a:2,X3,6i)} , (9) 

with the color factor Cp = 4/3, the number of colors Nc = 3, the symbol [dx] = dxidx2dx^ 
and the mass ratio r^r = rriQ/mB, ttLq being the chiral scale associated with the vr meson. 

The evolution factor Ei and hard function in Eq.Q can be found in Ref. [2^. In 
the derivation of spectator correction in the pQCD, we need to take the wave function of 
relevant mesons, we list the wave functions in appendix. 

For the decay, the CP asymmetry is time dependent, 



Acp{t) 



Sj,^^o sin(AMt) - cos(AMt) , (10) 



Where the mixing- induced asymmetry Sj/^^^o and direct CP asymmetry is defined as 

2 Im Aj/^^o 



S 
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1 + l^J/'iliwOl 

^J/^-^" - TVTx 12 ' \^^) 



where 



_ v;MJ/^p7r'\H,jf\-^) 

There are two ways to extract weak phase f3 through B^ ^ J /'^ tt" decay. The first way is 
through the dependence of the mixing-induced CP asymmetry on weak phase (3. The Sj/^^^o 
is not sensitive of input parameters, as shown in Fig. |H That means that the theoretical 
uncertainties of Sj/^.„o is quite small. If we measure the mixing- induced asymmetry Sj/^t,o, 
we can determine weak phase P through the dependence of Sj/^^^o on /5 as shown in Fig. [3] 
and Table [H 
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TABLE I: Determination of weak phase (3 through mixing-induced CP asymmetry Sj/^,^o 

Another way is to use the relation of the total asymmetry Acp with the weak phase /5. By 
integrating Acp{t)mth respect to the time variable t, we can get the total asymmetry Aqp, 

X „ 1 



A, 



CP 



1 +X2 



:S 



1 +X2 



:C 



(13) 



with X = Am/r ~ 0.723 for the B^-B^ mixing in the SM |21|. 

Like the mixing-induced asymmetry, the total asymmetry is also not sensitive to the 
input parameters, so we can determine the weak phase through the relation of the total CP 
asymmetry with weak phase (3 shown in Fig. [3l 

The numerical calculation needs some parameters and meson distribution amplitudes as 
input, we list them in the appendix. 

With the parameters and meson distribution amplitude in the appendix, we get the 
branching ratios oi B ^ J/"^ decays, AS'j/^^o and Cj/^^rO, 



[1.8a 



-0.182/ ,LN+0.0496/,,N+0.193^ p ^, +0.015 
-0.21 l'^"y'-0.02 \r)~0.171 

hO.866^, , +1. 173 \ +6.914 /\ +1.34 



0.014(/j/^)-o!o59 (A) -0.068 

(A)] X 10- 



Br{B° J/tpTT^) = 

Cj/^^o = [-9.936+--(^6)^--(7)^^;^^^(/.)li:?|(F0^^:^^(/5)] x 10' , 

= [2.84ttS^(^6)t°:3l(7);^ir(/^)lS:^^(i^i)tn^(/5)] x 10"^ (14) 

The main theoretical errors of the branching ratio are induced by the uncertainties below. 
The first error is from ub = 0.4 ± 0.04Ge\^, the second one is due to renormalization scale 
yU taken from mb/2 to mb, the third one is induced by 15% uncertainty of B ^ n form 
factor F^^"^, the fourth one arise from decay constant fj/^ = 0.405 ± 0.05GeV^, the fifth 
error is from CKM matrix parameter A = 0.2272 ±0.001, the sixth one is from CKM matrix 
parameter A = 0.818+°-°°^ 



-0.017- 



Compared with the experimental data 



2l| 



Br{B° J/tpn^) = (2.2 ± 0.4) x 10"^ 
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(15) 



our prediction of the branching ratio for 5 — J/\l/ vr^ is consistent with it. 

Unhke the branching ratio, ASj/^t^o and Cj/^t^o is not sensitive to CKM matrix parameter 
A or A, because these parameter dependences cancel out. The independence of ASj/^^o and 
Cj/i/)7r" on some CKM parameters is shown in Fig. |ll^a),(b),and Fig. [ni(a),(b). 

To find new physics and to extract the weak phase /3, we need rehable evaluation for the 
direct CP asymmetry Cj/^^o and AS'j/^^o, so we now consider the dependence of the direct 
CP asymmetry Cj/^t^o and ASj/^Trowith all parameters of input. The main uncertainties 
of Cj/^TjO and ASj/^t^o are induced by uncertainties of shape parameter ub, CKM matrix 
phase7, renormalization scale ^, B ^ n form factor Ff^'" and the weak phase p. The 
uncertainties of AS'j/^^o and Cj/^^^o are shown in Fig. |ll^c)-(f ) and Fig. [ni(c)-(f). 

Comparing with the result in Ref. 



Cj/^^o = 0.09 ±0.19 (16) 
Sj/^^o = -0.47 ±0.30 (17) 

our results of AS'j/^Troand Cj/^^o has much less theoretical uncertainties. So we conclude 
that if the measured deviation ASj/^t^o of the mixing-induced asymmetry is at 1% or the 
direct asymmetry Cj/^^o is at the level of percentage then we can say that there should be 
new physics . We are expecting precise measurement to the CP asymmetry of J / tpn^ 

in the near future. 



APPENDIX A: INPUT PARAMETERS AND WAVE FUNCTIONS 



We use the following input parameters in the numerical calculations 



Agr^) = 250MeV, U = ISOMeV, fe = 190MeV, 

iVLo 



rrin 



1.4GeV, Mb = 5.2792GeV, tbo = 1.53 x lO-^^s 



(Al) 



For the CKM matrix elements, we adopt the wolfenstein parametrization for the CKM 
matrix up to C(A^)|2l|, 

/ 



Vr 



CKM 



1 - f A AX^ip - 17]) \ 



-A 

\AX^{1- p-ir]) -AX^ 



1 - f AX' 



(A2) 
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with the parameters A = 0.2272, A = 0.818, p = 0.221 and r] = 0.340. 
For the B meson distribution amphtude, we adopt the model 2^ 



(A3) 



where ujb is a free parameter and we take Ub = 0.4 ± 0.05 GeV in numerical calculations, 
and Nb = 91.745 is the normalization factor for Uf, = 0.4. 

The J/ip meson asymptotic distribution amplitudes are given by 25j 



^^(x) = ^^(x) = 9.51 



fj/ip 
2v^ 



:X(1 — X) 



x{l — x) 



-\ 0.7 



x(l — X 



l-2.8x(l-x) 

n 0.7 



1 - 2.8x(l -x) 



^^(x) 



1.67 



2V2iVc 



l + (2x-l^ 



x(l — x) 



1 - 2.8x(l- x) 



0.7 



(A4) 



For the light meson wave function, we neglect the b dependant part, which is not impor- 
tant in numerical analysis. We choose the wave function of vr meson 

3 

u 



2^6 

fn 



2V6 



Ux{l -x)[l + 0.44Cf ^(2x - 1) + 0.25Cf ^(2x - 1^ 
"l + 0.43C2/^(2a; - 1) + 0.09C]/^(2x - 1^ 
1 - 2x) \l + 0.55(10a;2 - lOx + 1)" . 



(A5) 
(A6) 
(A7) 



The Gegenbauer polynomials are defined by 



Cl^\t) = U3t^ - 1), C'/\t) = i(35t4 - 30t^ + 3), 



C^/\t) = |(5t2 - 1), cr{t) = f{21t^ - lAe + 1). 
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FIG. 2: The dependence of the mixing-induced asymmetry Sjj^^o for J/'f Tr*^ on the 

weak phase [3 in diagram (a). The dependence of the deviation AS'j/^^o of the mixing- induced 
asymmetry from sin(— 2/3) on the weak phase (3 in diagram (6) 
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FIG. 3: The dependence of the the mixing- induced asymmetry Sj/^^o for J /'^ tt^ on the 

weak phase (3 in diagram (a) can be used to extract the weak phase (3 . The dependence of total 
CP asymmetry Acp on the weak phase (3 in diagram(6) can be used to extract the weak phase (3 
also. 
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FIG. 4: The uncertainties of A5j/^^o of the mixing-induced asymmetry from sin(— 2/3) are induced 
by that of renormahzation scale ^ in (c) , that of B ^ tt form factor in (d) , that of the weak phase 
7 in (e) and that of sin(2/3) in (/) . 
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FIG. 5: The uncertainties of the direct CP asymmetry Cji^^o are induced by that of renormahza- 
tion scale /i in (c) , that of i? — > vr form factor in (d), that of the weak phase 7 in (e) and that of 
sin(2/3) in (/) . 
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